Experimental estimation of the dimension of classical and quantum systems 
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An overwhelming majority of experiments in classical and quantum physics make a priori as- 
sumptions about the dimension of the system under consideration. However, would it be possible 
to assess the dimension of a completely unknown system only from the results of measurements 
performed on it, without any extra assumption? The concept of a dimension witness [1-5] answers 
this question, as it allows one to bound the dimension of an unknown classical or quantum system 
in a device- independent manner, that is, only from the statistics of measurements performed on it. 
Here, we report on the experimental demonstration of dimension witnesses in a prepare and measure 
scenario . We use pairs of photons entangled in both polarization and orbital angular momentum 
S to generate ensembles of classical and quantum states of dimensions up to 4. We then use 
a dimension witness to certify their dimensionality as well as their quantum nature. Our results 
open new avenues for the device-independent estimation of unknown quantum systems [sl-ITlll and 
for applications in quantum information science [l^, [13 • 
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Any physical theory aims at providing an explanation 
for the results of measurements performed on a system 
under different conditions. Given a system, one makes 
some general and plausible assumptions about its behav- 
ior, builds a theory compatible with these assumptions 
and tests its predictive power under different experimen- 
tal configurations. The dimension of the system, that is, 
the number of relevant and independent degrees of free- 
dom needed to describe it, is in general an a priori param- 
eter and represents one of these initial assumptions. In 
general, the failure of a theoretical model in predicting 
experimental data does not necessarily imply that the 
assumption on the dimensionality is incorrect, as there 
might exist a different model in a space of the same di- 
mension able to reproduce the observed data. 

A natural question is whether this approach can be 
reversed and whether the dimension of an unknown sys- 
tem, classical or quantum, can be estimated experimen- 
tally. Clearly, the best one can hope for is to provide 
lower bounds on this unknown dimension. Indeed, one 
can never exclude that more degrees of freedom are neces- 
sary to describe the system in more complex experimen- 
tal arrangements. The goal, then, is to obtain a lower 
bound on the dimension of the unknown system from 
the observed measurement data without making any as- 
sumption about the detailed functioning of the devices 
used in the experiment. Besides its fundamental inter- 
est, estimating the dimension of an unknown quantum 
system is also relevant from the perspective of quantum 
information, where the Hilbert space dimension is con- 
sidered as a resource (see, e.g., (l^). For instance, the 
optimal quantum realization of some protocols, e.g. bit 
commitment, requires systems of a minimal dimension 
[lij . Moreover, the dimension of quantum systems plays 
a crucial role in security proofs of standard quantum key 



distribution schemes that become insecure if the dimen- 
sion is higher than assumed [16]. 

The concept of a dimension witness allows one to es- 
tablish lower bounds on the dimension of an unknown 
system in a device-independent way, that is, only from 
the collected measurement statistics. It was first in- 
troduced for quantum systems in connection with Bell 
inequalities in and further developed in [H, [l7l - [22| . 
Other techniques to estimate the dimension have been 
developed in scenarios involving random access codes ^ , 
or the time evolution of a quantum observable 

More recently, a general framework for the study of 
this question has been proposed in [5|. In this approach, 
dimension witnesses are defined in a prepare and measure 
scenario where an unknown system is subject to different 
preparations and measurements. One of the advantages 
of this approach is its simplicity from an experimental 
viewpoint when compared to previous proposals. The 
considered scenario consists of two devices (see Fig. [1]), 
the state preparator and the measurement device. These 
devices are seen as black boxes, as no assumptions are 
made on their functioning. The state preparator pre- 
pares upon request a state. The box features N buttons 
which label the prepared state; when pressing button x, 
the box emits a state where x G {1, ...,n}. The pre- 
pared state is then sent to the measurement device. The 
measurement box performs a measurement y G {!,..., m} 
on the state, delivering outcome b G {!,..., /c}. The ex- 
periment is thus described by the probability distribution 
P{b\x^y)^ giving the probability of obtaining outcome b 
when measurement y is performed on the prepared state 
px' The goal is to estimate the minimal dimension that 
the ensemble {px} must have to be able to describe the 
observed statistics. Moreover, for a fixed dimension, we 
also aim at distinguishing sets of probabilities P{b\x^y) 
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FIG. 1: Sketch of a device- independent test of the dimen- 
sion of an ensemble of quantum states. Our scenario features 
two black boxes: a state preparator, which generates one of n 
possible states on demand, and a measurement device which 
performs one of m possible measurements on the state gener- 
ated. 

that can be obtained from ensembles of quantum states, 
but not from ensembles of classical states. This allows 
one to guarantee the quantum nature of an ensemble of 
states under the assumption that the dimensionality is 
bounded. 

Formally, a probability distribution P{b\x^y) admits a 
(i- dimensional representation if it can be written in the 
form 

P{b\x,y)=tr{p^M^), (1) 

for some state px and operator acting on C^. We 
then say that P(6|x, y) has a classical d-dimensional rep- 
resentation if any state of the ensemble {px} is a classical 
state of dimension i.e., a probability distribution over 
classical dits (or equivalently, all the states of the ensem- 
ble commute). 

A dimension witness for classical (quantum) systems 
of dimension d is defined by a linear combination of the 
observed probabilities P{b\x^y)^ defined by a tensor of 
real coefficients Di)^x,y^ such that 

^Db,x,yP{b\x,y)<Wd (2) 

b,x,y 

for all probabilities with a d-dimensional representation, 
while this bound can be violated by a set of probabili- 
ties whose representation has a dimension strictly larger 
than d. A family of dimension witnesses In was intro- 
duced recently [5| in a scenario consisting of TV possible 
preparations and — 1 measurements with only two pos- 
sible outcomes, labeled by b = ±1. Here we shall focus 
on one of these witnesses, namely 

I4 = Ell + E12 + Ei3 

+ E21 + E22 — E23 

+ E31 — 

- ^41, (3) 

where Exy = P{b = -\-l\x^y) — P{b = -\-l\x^y). Witness 
I4 features four preparations and three measurements. It 





C2 (bit) Q2 (qubit) C3 (trit) Q3 (qutrit) C4 (quart) 


h 


5 6 7 7.97 9 



TABLE I: Classical and quantum bounds for the dimension 
witnesses I4. The maximal possible value of I4 for classical 
systems of dimension d is denoted Cd, while the maximum 
possible value for quantum systems of dimension d is denoted 
Qd' The witness can be used to discriminate ensembles of 
classical and quantum states of dimension up to d = 4. Note 
that for a given value of the dimension of the system, the 
witness can be used to certify its quantum nature. 



can distinguish ensembles of classical and quantum states 
of dimensions up to d = 4. All the relevant bounds are 
summarized in Table HI 

In order to test this witness experimentally, we must 
generate classical and quantum states of dimension 2 
(bits and qubits, respectively), classical and quantum 
states of dimension 3 (trits and qutrits), and classical 
states of dimension 4 (quarts). To do so we exploit the 
angular momentum of photons 0, . The total angular 
momentum contains a spin contribution associated with 
the polarization, and an orbital contribution associated 
with the spatial shape of the light intensity and its phase. 
Within the paraxial regime, both contributions can be 
measured and manipulated independently. The polar- 
ization of photons can be described in a 2-dimensional 
Hilbert space, as any polarization state can be described 
by a weighted superposition of two orthogonal polariza- 
tion states (e.g., horizontal and vertical). The spatial 
degree of freedom of light lives in an infinite-dimensional 
Hilbert space, as any spatial waveform can be described 
in terms of a weighted superposition of modes that span 
an infinite-dimensional basis [23[. Paraxial Laguerre- 
Gaussian (LG) modes constitute one of these bases. LG 
beams are characterized by two integer indices, p and m. 
The index m determines the azimuthal phase dependence 
of the mode, which is of the form ~ exp {imcp)^ where cp 
is the azimuthal angle in the cylindrical coordinates. LG 
laser beams carry a well-defined orbital angular momen- 
tum (0AM) of mh per photon that is associated with 
their spiral wavefronts |24| . 

Our state preparator employs a source based on spon- 
taneous parametric down-conversion (SPDC) that gen- 
erates pairs of photons (signal and idler) entangled in 
both polarization and 0AM. By performing a projective 
measurement on one photon of the pair (idler), we pre- 
pare its twin photon (signal) in a well-defined state. A 
given measurement outcome on the first photon (idler) 
corresponds to a given preparation for the second pho- 
ton (signal), which represents the mediating particle of 
Fig. 1. 

We use the state preparator to generate ensembles of 
states of dimension up to 4: two degrees of freedom 
are encoded in their polarization (horizontal or verti- 
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FIG. 2: Experimental setup of the dimension witness. The 
upper part depicts the individual states \(l)x) {x = 1...4) pre- 
pared by the state preparator. The lower part illustrates the 
measurements performed by the measurement device. SF - 
spatial filter; L - lens; M - mirror; XL - nonlinear crystal; DL 

- delay line; F - 5-nm bandpass filter; BS - beamsplitter; HWP 

- half-wave plate; P - polarizer; SLM - spatial light modula- 
tor; PBS - polarizing beamsplitter; FL - Fourier-transform 
lens; SMF - single-mode fiber; D - single-photon counting de- 
tector; & - coincidence electronics. 



cal) and two degrees of freedom are encoded in their 
0AM (m = +1 or m = —1). For the case of qubit 
preparations, the four states l^a^) {x = 1...4) are en- 
conded in the polarization only and have the same 0AM 
of m = -hi (see the upper left part of Fig. [2j). An 
analogous procedure is used to prepare a qutrit state, 
however, now state 1^3) carries 0AM of m = —1, thus 
adding one dimension. The quart state is prepared by 
flipping the 0AM of state {(^a). In this way we gen- 
erate a 4-dimensional Hilbert, spanned by the orthogo- 
nal vectors {\H, +1), \H, -1), \V, +1), and|V, -1)}, where 
±1) and | V, ±1) stand for a horizontally or vertically 
polarized photon with 0AM equal to m = ±1, respec- 
tively. 

To implement continuous transition from quantum to 
classical states, a polarization-dependent temporal delay 
between the two photons of an entangled pair is intro- 
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FIG. 3: a) Computer generated hologram used to detect a 
LG beam with m = +1. b) CCD camera photograph of the 
generated LG mode when the hologram is illuminated by a 
Gaussian beam. 



duced. If the temporal delay between the photons ex- 
ceeds their correlation time, the coherence is lost, i.e., 
the off-diagonal terms of the state in the ensemble go to 
zero (see Appendix). 

Finally, the measurement device simply performs mea- 
surements on the second (signal) photon. 

The experimental setup is shown in Fig. [2l The second 
harmonic (Inspire Blue, Spectra Physics/Radiantis) at a 
wavelength of 405 nm of a Ti:saphire laser in the picosec- 
ond regime (Mira, Coherent) is shaped by a spatial filter 
and focused into a 1.5-mm thick crystal of beta-barium 
borate (BBO), where SPDC takes place. The nonlinear 
crystal is cut for collinear type-II down-conversion so that 
the generated photons have orthogonal polarizations. Be- 
fore splitting the photon pairs towards the state prepara- 
tor and measurement device, the polarization-dependent 
temporal delay r is introduced. The delay line (DL) con- 
sists of two quartz prisms whose mutual position deter- 
mines the difference between the propagation times of 
photons with different polarizations. 

The measurement setup consists of a half-wave plate 
(HWP), polarizer (P), spatial-light modulator (SLM) 
and a Fourier-transform lens (FL). The half- wave plate 
and polarizer project the incoming idler photon into the 
desired polarization state. The desired 0AM state is se- 
lected by the spatial light modulator. If the nonlinear 
crystal is pumped by the fundamental LG beam LGqo, 
about 15% of the entangled photons are generated in the 
LG mode with m = ±1 125|. The 0AM is conserved 



28| (see Appendix for more details), i.e., if the sig- 
nal photon has m = +1, the idler must have m = — 1 
and vice versa. The SLM encodes computer-generated 
holograms that transform the m = +1 state or m = — 1 
state into the fundamental LG state LGqo [6] that is cou- 
pled into a single- mode fiber (SMF). In this way, only 
a photon having a correct combination of its 0AM with 
a corresponding hologram placed at the SLM is success- 
fully coupled in the fiber and results in a subsequent click 
at the detector. Figure [3] a) shows a computer-generated 
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Delay t (fs) Hologram position (arb.u.) 

FIG. 4: a) The Hong-Ou-Mandel dip shows the transition 
from a pure entangled state to a mixed state. The pure state 
is used to prepare a quantum superposition of \H) and \V), 
whereas the mixed state is used to prepare statistical mixtures 
of \H) and \V). b) Two-dimensional scan of the hologram 
position in the signal beam. Signal hologram is set to m = 
+ 1 and idler hologram is set to m = 0. Disappearance of 
coincidences asserts the correct position of the hologram was 
found in order to detect signal photons with 0AM m = +1. 

hologram to measure an m = +1 state and Fig. [3] b) 
shows a picture of the LG mode taken by a GGD camera 
when this hologram is illuminated by a Gaussian beam. 

When the polarizers in both arms are set to +45 de- 
grees, a Hong-Ou-Mandel dip [29] of visibility above 95% 
is measured (see Fig. [Ha)). The center of the dip corre- 
sponds to the pure state generation when the qubit and 
qutrit states are prepared depending on their 0AM. On 
the right-hand side, on the shoulder of the dip, the su- 
perposition is lost and the qubit and qutrit are converted 
to a classical bit and trit, respectively. The purity of the 
state is given by Eq. [8] that can be found in Appendix 
with more details. 

To find the correct position of the hologram on the 
SLM, the center of the hologram must exactly coincide 
with the center of the down-converted beam. To overlap 
them perfectly, first the signal hologram is set to detect 
photons with m = +1 and the idler hologram is set to de- 
tect m = 0. Then signal holograms with different central 
positions of the vortex are scanned across the SLM. Since 
the down-converted photons are generated with m = ±1, 
a minimum of coincident counts between m = +1 and 
m = tells us the correct position of the signal hologram 
(see Fig. IHb)). Now with the signal hologram fixed, the 
idler hologram set to m = —1 is scanned to find a max- 
imum number of coincidences. With both holograms in 
the correct place, the dimension witness is measured. 

First, the qubit ensemble is generated. We sequentially 
prepare the four qubit states l^^^) of Fig. [2] and perform 
the three desired measurements. The first measurement 
y = 1 assigns dichotomic measurement results of 6 = +1 
and 6 = — 1 to horizontally and vertically polarized pho- 
tons, respectively. The second measurement ^ = 2 as- 
signs 6 = +1 and 6 = — 1 to 0AM values of m = +1 and 
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FIG. 5: Dimension witness I4 for qubit (blue triangles), qutrit 
(red circles) and quart (green diamonds) as a function of tem- 
poral delay r. For delays r > 255 fs, coherence is lost and 
quantum superpositions turn into statistical mixtures, i.e., 
classical states. The maximum violations for qubit, qutrit 
and quart are 5.66 di 0.15, 7.57 ± 0.13, and 8.57 =b 0.06, re- 
spectively, which are close to the corresponding theoretical 
bounds, given in Table H] which are represented here by the 
horizontal lines. 



m = — 1, respectively. And the third measurement y = 3 
assigns 6 = +1 and 6 = — 1 to photons polarized at +45^ 
and —45^. The maximum theoretical value of the witness 
I4 of Eq.[3]for this ensemble of states and this set of mea- 
surements is 5.83 (see Appendix). The experimentally 
measured witness reaches I 4 = 5.66 ± 0.15. This clearly 
demonstrates the quantum nature of our 2-dimensional 
system, since classical bits always satisfy I4 < 5. 

Now a variable temporal delay r is introduced between 
the horizontal and vertical polarizations, and the value 
of the witness decreases below 5, leveling off at r 255 
fs that corresponds to the end of the Hong-Ou-Mandel 
dip, as expected (see blue triangles of Fig. [5j). Thus our 
qubit is now converted into a classical bit. 

If the 0AM of state \(j)s) is flipped, a qutrit is prepared 
and the dimension witness is measured. For a zero delay 
r, the value of the witness reaches I4 = 7.57 ± 0.13, cer- 
tifying the presence of a (at least) 3-dimensional quan- 
tum system. Here the maximum theoretical value for 
this ensemble of states and this set of measurements is 
5 + 2^2 - 7.83 (see Appendix). 

When the delay line is scanned, the values of the wit- 
ness drop below 7, but in a certain range of r they remain 
above the qubit bound of 6, testifying that at least three 
dimensions are present (see red circles of Fig. [5]). In the 
classical limit of large delays, the witness values are still 
larger then the bound of /4 = 5 for bits. 

Finally, we prepare classical 4-dimensional systems, i.e. 
quarts. Now the first measurement y = I assigns the di- 
chotomic value of 6 = +1 to horizontally polarized pho- 
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tons of m = ±1 and to vertically polarized photons of 
m = +1, and b = —1 corresponds to vertically polarized 
photons with m = — 1; the second measurement y = 2 
assigns 6 = +1 and 6 = — 1 to horizontally and vertically 
polarized photons, respectively; and the third measure- 
ment y = 3 assigns 6 = +1 and b = —1 to 0AM values 
of m = +1 and m = — 1, respectively. All the correla- 
tion coefficients should add up and should be equal to 1. 
The obtained value is /4 = 8.57 ±0.06, which violates the 
qutrit threshold by more than 10 standard deviations. In 
this case, the values of the witness remain constant and 
are independent of the temporal delay r. Indeed this is 
because the state is classical (a statistical mixture of or- 
thogonal quantum states) and no superposition is present 
(see green diamonds of Fig. [5j). 

To conclude, we have demonstrated how the dimension 
of classical and quantum systems can be bounded only 
from experimental measurement results without any ex- 
tra assumption on the devices used in the experiment. 
Dimension witnesses represent an example of a device- 
independent estimation technique, in which relevant in- 
formation about an unknown system is obtained solely 
from the measurement data. Device-independent tech- 
niques (see also Refs. [sl-llll) provide an alternative and 
robust approach to the estimation of quantum systems, 
especially when compared to most of the existing esti- 
mation techniques, such as quantum tomography, which 
crucially rely on some assumptions about the measured 
system, e.g., its Hilbert space dimension, that may be 
questionable in complex setups. Our work demonstrates 
how the device-independent approach can be employed 
to experimentally estimate the dimension of an unknown 
system. 



APPENDIX 

In order to generate photons in the multidimensional 
Hilbert space defined by the angular momentum of light, 
and at the same time to be able to modify the dimension 
and the properties of the corresponding quantum states, 
we generate photon pairs entangled in the polarization 
and spatial degrees of freedom [7]. 

If the beam sizes of the pump beam and those of 
the collection systems of the down-converted photons are 
large enough, we can neglect the effects of the Foynting- 
vector walk-off that takes place in a birefringent medium 
0. In this case, the generated paired photons are found 
to be embedded into LG modes with opposite winding 



numbers m. 



-rrij, where m. jh refers to the 0AM 



of the signal and idler photons, respectively. After the 
down-conversion stage, a beam splitter is used to sepa- 
rate the signal and idler photons. If the optical system 
selects only the modes with rris^i = ±1, the quantum 
state of the photon pairs can be written as 



= ^ [ dQ 

X {f{n)[\n,H,^i),\n,v,-i)2 

+ |l],i7,-l)l|l],F,+l)2] 

+ f{-n)[\n,v,^i)i\n,H,-i)2 

+ |1],F,-I)l|l],i7,+1)2]}. 



(4) 



The signal photon has frequency uJs = oJq -\- and the 
idler photon has frequency uji = ujq — where Uq is the 
central frequency of the down- converted photons and Q is 
the frequency deviation from the central frequency. The 
function / reads 



m) 



_sincf^jexp(i^j, (5) 



where L is the crystal length, and D = 1/vs — 1/vi is the 
difference of inverse group velocities (vs^i) of the signal 
and idler photons. 

The signal and idler photons are entangled in the po- 
larization and spatial degrees of freedom. On the other 
hand, the frequency properties of entangled two-photon 
states cannot be neglected even when the entanglement 
resides in the other degrees of freedom. In most ex- 
periments, the entangled states take advantage of only 
a portion of the total two-photon quantum state and 
the generation of high-quality entanglement requires the 
suppression of any frequency "which-path" information 
that would otherwise degrade the degree of entanglement. 
Notwithstanding, here we use the presence of frequency 
"which-path" information to tailor the quantum nature 
of ensembles of signal photons. 

We introduce a temporal delay r between the orthog- 
onally polarized signal and idler photons before they are 
separated at the beam splitter. In the idler path, a po- 
larizer rotated by an angle cp^ is placed, followed by a 
spatial light modulator which in combination with a sin- 
gle mode fiber selects idler photons in a LG mode with 
a winding number m = ±1. Thanks to the angular mo- 
mentum correlations between signal and idler photons, 
the detection of an idler photon with horizontal or verti- 
cal polarization {cp = 0^) and with a spatial shape given 
by a LG mode with index m = ±1, projects the signal 
photon into the orthogonal polarization and a LG mode 
with index m = ^1. When the polarizer projects the 
idler photon in a weighted combination of the two or- 
thogonal polarizations {cp 7^ 0^), the resulting quantum 
state of the signal photon depends on the delay r. 

The Hilbert space of signal photons, defined by 
the angular momentum (i.e., polarization and or- 
bital angular momentum), is spanned by vectors 
{\H, +1), \H, -1), \V, +1), and|F, -1)}. After tracing out 
the frequency degree of freedom, the quantum state of 
signal photons is given by the density matrix 
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a cos^ (p 

aY{r)sm2(p/2 

{l-a)Y{r)sm2ip/2 



a-f{r) sin 2if/ 2 \ 

(1 - a)-f{r) sin 2(p/ 2 

a sin^ if 



(6) 



where a = 1 if the idler photons is projected into a mode 
with m = — 1, and a = if the idler photon is projected 
into a mode with m = 1. 

The parameter 7(r), which fulfills the condition 
|7(t)| < 1, writes 
DL 
~2^ 



7(r) 



/ 



dVtsinc^ exp {iVt (2r 



DL)} 



tri 



DL/2^ 



max(l 



DL/2 

where tri(t) is the triangle function tri(t) 
|t|,0). 

The purity of the state given by Eq. ([6]), provided a 
1, writes 

l7(r)P 



(7) 







1 



1 



sin 2Lp. 



(8) 



The parameter 7(r) defines the "quantumness" of the 
ensemble of states. A d-dimensional ensemble of states 
{px} = ^ ' ' ' d)^ where all the states satisfy the condi- 
tion 7a: (r) = 0, constitutes an ensemble of classical states, 
while if any of the states has 7(r) ^ 0, the ensemble is 
of quantum nature. 

The quart is generated by the four combinations of 
a = 0, 1 and Lp = 0,90^. A three-dimensional quantum 
system (qutrit) is generated by states : a = l^(p = 
22.5^ \^2) : a = l,if = -22.5^ \h) : a = 0, = 0^ 
and 1^4) : a = l^cp = 90^. The delay line is set to 
r = DL/2 so that |7^| = 1. When \-fx{r)\ 0, the 
qutrit is converted into a classical trit. The theoretical 
value of I4 for this ensemble of states and this set of 
measurements is 



/4 = 5 + 2 cos 2(p + 27(r) sin 2(p 



(9) 

that attains its maximum value of 5-\-2y^l + |7(t)P when 
tan2(/:) = 7(r). For 7 = 1, we obtain /f = 7.83 > 7. 
For 7 = 0, the maximum value a classical trit can reach 
is 7. 

A two-dimensional quantum system (qubit) is gener- 
ated by the same set of states as in the case of qutrit, 
except that state |c) now has a = 1 and cp = 0'^. A clas- 
sical two-dimensional system (bit) is obtained for 7 = 0. 
The value of is 

/4 = 3 + 2 cos 2(/:) + 27(r) sin 2cp (10) 

with its maximum of 3 + 2^/1 + |7(t)P. For 7 = 1, 
we obtain /^^^ = 5.83 > 5 and for 7 = we obtain 

jmax _ ^ 
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